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ABSTRACT 

We investigate the main physical properties of low-metalhcity Asymptotic Giant 
Branch stars, with the aim of quantifying the uncertainties that presently affect the 
predicted chemical yields of these stars, associated to mass loss and description of 
molecular opacities. We find that above a threshold mass, M ~ 3.5 Mq for Z = 0.001, 
the results are little dependent on the opacity treatment, as long as hot-bottom burn- 
ing prevents the surface C/0 ratio from exceeding unity; the yields of these massive 
AGB stars are expected to be mostly determined by the efficiency of convection, with 
a relatively mild dependence on the mass-loss description. A much higher degree of 
uncertainty is associated to the yields of less massive models, which critically depend 
on the adopted molecular opacities. An interval of masses exists, say 2.0 — 3.0 M©, (the 
exact range depends on mass loss) , in which HBB may be even extinguished following 
the cooling produced by the opacity of C-bearing molecules. The yields of these stars 
are the most uncertain, the variation range being the largest (up to ^ 2 dex) for the 
nitrogen and sodium yields. For very low-mass models, not experiencing hot-bottom 
burning (Af ^ 1.5 Mq), the description of mass loss and the treatment of the con- 
vective boundaries are crucial for the occurrence of the third dredge-up, with sizable 
consequences on the CNO yields. 
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1 INTRODUCTION 

The Asymptotic Giant Branch (AGB) is a common phase 
of the ev olution of stars with mass in the range l-SM© 
(|Busso et al. 1999i : lHerwieil2005l ). After core- helium burning, 
they experience a series of thermal pulses (TPs), triggered 
by the ignition of helium in a thin layer below the CNO 
burning shell, under conditions of thermal instability. The 
strong mass loss that they experience favours the ejection 
of the whole external mantle, with the formation of a CO 
white dwarf. 

Historically, the interest towards AGBs s tarted with th e 
discovery of s-process enriched AGB stars (|Merril]||l952f ). 
and with their identification as sites of nucleosynthesis and 
s-process (|Burbidge et al.|[l957l '). and was futher stimulated 
by the discovery of carbon and lithium-rich stars in the 
Galajcy and the Magellanic Clouds (MCs) (Blan co et al.l 
ll980l : lAbiaet al.| [l991: Smith & Lambert 1989, 199^): they 
were correctly interpreted as thermally pulsating stars 
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ISackmann fc Boothrovdl Il992l : llben fc Renzinilll982l '). Ex- 
tended investigations of the physical properties of AGB 
stars showed the existence of two mechanisms able to al- 
ter the surface chemistry: the Third Dredge-Up (TDU), 
i.e. the inwards penetration (in mass) of the external en- 
velope after ea ch TP into regions earlier involved in 3a 
nucle osynthesis (|lbenlll975l : lFuiimoto. Nomoto. &: Sugimotd 
Il976l ): and Hot Bottom Burning (HBB), when the base 
of the convective envelope becomes sufficiently hot (^ 
40 X lO^K) to produce pr oton-cap ture nuc leosynthesis 
(Sc alo. Despain. fc Ulrichl [l975: Rcnz ini fc Volil il981'). and 
enough energy to cause the brea kdown of the core mass- 
luminosity relation on the AGB ([Blocker fc Schoenberneij 
ll99ll : lBoothrovd fc Sackmannlll99a V 

In more recent years, we have been gathering the com- 
pelling evidence that TP-AGB stars play a crucial role in 
many properties of their host systems, mainly owing to their 
intrinsic brightness and distinctive spectral features. For in- 
stance, the TP-AGB contribution to the total luminosity 
of single-burst stellar population s reaches a maximum of 
about 40% at ages from 1 to 3 Gyr (jFrogel. Mould, fc Blanco! 
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Il990r ). and accounts for most of the bright-infrared ob- 
jects in resolved galaxies, as clearly demonstrated by 
DENIS, 2MASS SAGE, S^MC, and AKARI IRC data 
JCioni et al.lll999l: 



l2006l : lBolattoet al 



Nikolaev fc Weinberd [20001 : iBlum et alj 



20071 : llta et al.ll2008l ') for the Magellanic 



Clouds (MC). Observational data of AGB stars is now avail- 
able also for other galaxies of the Local Group - su ch as Leo I 
iJHeld et al."2010'). Leo II (GuUicuszik et al. 2008), Sagittar- 
ius (iGuUieuszik et al. 2007). P hoenix (Menzies et al. 2008 ). 



Fornax (jWhitelock et al.ll2009l ) -, as well as for more distant 
dwarf galaxies (dalcanton 2008]j. 

Chemical pollution from massive AGB stars is cur- 
rently one of the most plausible hypothesis to account for 
the c hemical pattern s observed in Globular Cluster (GC) 
stars (|Carrettall2006r ) . w here they may explai n the existence 
of multiple populations (jVentura et al.ll200ll ). indicated by 
t he discovery of multiple main sequences in some GCs 
Piotto et ai]l2007l ). The pioneering model bv lD'Ercole et all 



20081 ) outlines the capability of massive AGB stars to pro- 



vide an efficient pollution of the interstellar medium, and to 
stimulate the formation of new stellar generations, with pre- 
dicted photometric and spectroscopic properties consistent 
with observations. 

Moving from zero to high redshift, the light contri- 
bution from AGB stars becomes significant in the opti- 
cal, where galaxies are dominated by intcrmcdiate-age stars 
jBressan. Chiosi. fc Fagotto 1994 : Bruzual fc Chariot 2003= : 
iMarastonI |2005| ). This fact is crucial: it has been recently 
pointed out that quantifying the weight of TP-AGB stars 
has a lar ge impact on the mass assembly in high redshift 
galaxies (jSantini et al.l bOOO'). 

The prominent role played by AGB stars in so many 
astrophysical contexts has increased the demand for de- 
tailed AGB modelling. Since full computations of the TP- 
AGB phase are extremely time-consuming, given the short 
time steps that must be used during each TP, many syn- 
thetic models have been developed and are presently used, 
in which, based on existing full models, the AGB properties 
are parametrized as a function of stellar mass and mctallic- 
ity (jy an den Hock & Grocncwcgcn 1997; Marigo ct al. 199_9i; 
llzzard et al.ll2004l : iMarigo fc Girardill2007i ). 

Unfortunately full AGB modelling is dependent on the 
assumptions on many physical ingredients that are not 
known from first principles, and must be described by means 
of semi-empirical calibrations. The treatment of convection 
is ascribed most of the large differences among AGB mod- 
els in the existing literature: the efficiency of the convective 
instability, and the assumptions concerning the extension 
of the mixed regions have a strong impact o n the physical 
and chemical evolution of massive AGB stars l|Herwieil200d . 
l2005l : IVentura fc D'Antona|[2005al ) . 

Recent investigations have also shown that the main 
evolutionary properties of AGB stars can be critically af- 
fected by the radiative low-T opacities adopted to model 
their outer mantles. The use of opacities that correctly ac- 
count for the drastic changes in the molecular chemistry of 
the gas, when the C/O ratio passes from below to above 
unity, is mandatory, otherwise the general expansion of the 
outer layers triggered by the opacity increase is missed 
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(jMarigol [2003). Since one of the most significant conse- 
quences of a correct computation of the molecular opacities 
is a sudden increase in mass loss, a careful study of the effects 
driven by molecular opacities on the evolutionary properties 
and stellar yields requires to investigate the impact of mass 
loss, and its uncertainties. 

In order to produce more reliable full AGB models, that 
may also contribute a useful input to synthetic AGB models, 
we have started a new project dedicated, in a first phase, to 
understand and quantify the main uncertainties associated 
to the various input physics. 

The properties of the most massive AGB stars, and the 
rel evant role played by convectio n modelling, was analysed 
by Ventura fc D'Antona l|2005al ). Moving to lower masses, 
IVentura fc Marigol (2009) have discussed the role played by 
molecular opacities. In this work we concentrate on how the 
opacity treatment is interfaced with the description of mass 
loss, and which results in the literature need a substantial 
revision. Future steps will be an extensive comparison of 
model predictions with observations, to provide the astro- 
physical community with more reliable yields. 

The paper is organized as follows. Section[2]summarizes 
the uncertainties affecting AGB models. The physical ingre- 
dients of the models presented in the present investigation 
are detailed in Sect. [3] Section|4]is dedicated to the analysis 
of the physical properties of the models, and how they de- 
pend on the choices concerning opacity and mass loss. The 
reliability of the yields, and the expected chemical patterns, 
are discussed in Sect. [5] A comparison with the results from 
other investigations on the same topic is presented in Sect.[Bl 



2 ABOUT THE UNCERTAINTIES 
AFFECTING AGB MODELLING 

The theoretical description of the AGB evolution is ex- 
tremely sensitive to the input physics used to calculate the 
models: the treatment of the convective instability, the de- 
scription of mass loss, the extent of the possible extra- mixing 
region beyond the formal convective/radiative boundary de- 
termined by the classic Schwartzschild criterion, the adopted 
low-temperature opacities, and the nuclear cross-sections of 
some p- and a — capture react ions. An exhaustive analysis 
can be found inlHerwid (|2005|). 

Ventura fc D'Antona (l2005al ). f ollowing the ana l ysis b y 



'Blocker fc Schoenberner' (1991) and lRenzini fc Volil (|l98lh . 
showed that convection modelling plays an important role in 
determining the evolution during the TP phase: the differ- 
ences in the convection modelling have been recognized as 
the main cause responsible for the d iscrepancies among the 
resul t s obtained by different groups (Denisscnko v fc Herwid 
l2003l : iFenner et a"l]r2004 : Karakas fc Lattanzio 20071 ). 

The assumption of some convective overshoolu from 
the bottom of the envelope is closely associated to the oc- 
currence of TDU after a TP that, in turn, is confirmed 
by the discovery of many carbon s tars in the Galaxy 
and in t he Magellanic Cloud s (see e .g. iBattinelli fc DemersI 
(|2004l) : iGroenewegen et al.1 (|2009l )). The analyses by 



■^ By convective overshoot we mean here the further distance 
travelled by convective eddies beyond the border where buoyancy 
vanishes, fixed by the Schwarzschild criteriuni 
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iGroenewegen fc de JoiiS (|l993l ') and iMarigo et al] (|l999l ). 
based on synthetic TP-AGB models, indicate that the extent 
of the TDU required to reproduce the luminosity function 
of the MCs should be larger than predicted by the stan- 
dard models, where TDU is found in the context of the 
Schwarzschild criterion. 

With respect to mass loss along the AGB, several pre- 
scriptions, based on either theoretical or empirical grounds, 
have been proposed in the literature. Among them we re- 
call: Reimers' classical law (1975) wa s a fit of measured 
mass- loss rates as a function of MR/L; IVassiliadis fc WoodI 
([1993) (hereinafter VW93) was based on the observed cor- 
relation between mass-lo s s rat es and periods of pulsat- 
ing AGB stars; iBlockeiJ (|l995l ) proposed a modi fication 
to the Reimers' law base d on fi tting models from iBowenI 
l|l98g) : iBowen fc Willsonl ([199l|) presented a prescription 
based on dynamical models of pulsating oxygen-rich atmo- 
spheres includi n g the formation of dust; basing on those 
models iBediinI ([1983) derived a mass-loss formalism as 
a function of basic stellar parameters, Af , R, and Tctt; 
IWachter et al.l (|2002l . 120081 1 proposed a fo rmula based on 
dust-d riven superwind models of C stars; Ivan Loon et al.l 
(|2005|) presented an empirical calibration based on obser- 
vations of d ust-enshrouded red superg iants and oxygen-rich 
AGB stars; ISchroder fc Cunt j ([20051) introduced a correc- 
tion to the Reimers' fo r mula based on some physical argu- 
ments; IStraniero et al.l (|2006l ) (hereinafter S06) formulated 
a revised calibration of the mass loss-period relation for pul- 
sating AGB stars. 

In AGB models, an additional source of uncertainty 
is related to the adopted low-T radiative opacities, which 
quantify the removal of radiation energy as photo ns pass 
throug h the H-rich mant le. The early predictio ns bv lMarigol 
(|2002h were confirmed bv lCristallo et al.l (|2008l ). who showed 
how the physical and chemical evolution of low-mass, low- 
Z stars are affected by the adoption of the correct opacity 
treatment. 



3 PHYSICAL AND CHEMICAL INPUTS 

The models presented in this paper were calc ulated by 
means of the ATON code for stellar evolution (|Mazzitellil 
19791). a full desc ription of which can be found in 



Ventura et all (|l998l ). 



3.1 Convective regions 

Convection was modelled according to the Full Spec- 
trum of Turbulenc e (her einafter FST) prescription by 
ICanuto fc Mazzitellil (|l99lh . Mixing of chemicals and nu- 
clear burning were coupled by mea ns of a d iffusive approach, 
following the scheme presented in ICloutm an & EoU (1976); 
accordingly, convective overshoot was modelled by an ex- 
ponential decay of convective velocities beyond the formal 
borders, with an e- folding decay of I — C,Hp. During the 
two main phases of core nuclear burning and in occurrence 
of the second dredge-up we used C, — 0.02, in agreement with 
the calibration based on the width of the m ain sequences of 
open clusters given in lVentura et al.l (|l998l ). During the TP- 
AGB phase, given the uncertainties outlined in the previous 
section, no overshoot was considered. 



Figure 1. Evolution of the mass-loss rate {M in Mq/jt) as a 
function of the surface luminosity for models with initial masses 
of SMq calculated with different prescriptions for mass loss and 
molecular opacity (see Table [T] for the meaning of the various 
symbols). Each point in the plot marks the quiescent stage of 
the luminosity maximum preceding the occurrence of a thermal 
pulse. The solid lines connect the four evolutionary sequences at 
the stage when the total mass of the star has been reduced to 2.9, 
2.7, 2.5 and 2.0 Mq. 



3.2 Mass loss 

To investigate the sensitivity of the results obtained on 
the mass loss prescription, our choice was to compare the 
findings obtain ed w ith two different presc riptions, namely 
iBlockerl (liggsl ) and IStraniero erahl (|2006l '). In the former 
case it is ass umed, based on hydrodynamical simulations 
(|BowerJll988h . a steep increase of mass loss with luminos- 
ity as the star enters the AGB phase. This is simulated by 
multiplying the canonical Reimers' formula by a luminos- 
ity power (L^'^). The free parameter entering the Reimer's 
prescription is rjR = 0.02, according to the calibration of 
the luminosity function of l ithium-rich sta r s obse rved in the 
Magell anic Clouds give n in Ventura et al.l ([20001). 

In IStraniero et al.l (|2006r ) the mass-loss rate is made 
vary with the pulsation period in the fundamental mode: 
compared to VW93 the empirical revision by S06 (see their 
figure 5) , based on a compilation of more recent data, pre- 
dicts larger mass-loss rates for pulsation periods in the range 
100 — 300 days, whilst lower rates are associated to periods 
larger than 500 days. 



3.3 Equation of state & opacities 

For the equation of state (EOS) we adopted the latest ver- 
sion of the OPAL EOS (200 5) where available, superseded 
by the ISaumon et al.l ([1995|) EOS in the partial ionization 
regime. The EOS is extended to the high-density, high- 
temperature regime a c cordin g to the description given in 
IStoltzmann fc Blocker! (|2000l '). 

The radiative opacities were calculated follow i ng th e 
OPAL prescription, according to llglesias fc Rogers! ([1999). 
The conductive opacities were taken from Poteckin (2006 jj. 

As to the low-temperature opacities, i.e. 1 500 K ^ T ^ 
10 OOP K, we used t he sa me large data set of tables as in 



IVentura fc Marigol (|2009|), which were computed with the 
iESOPUS tool (Marigo & Aringer 2009 j3 to consistently fob 
low the significant changes in the CNO surface abundances 
caused by the TDU and HBB during the TP-AGB evolution. 
The opacity set was designed to account for the com- 
plex interplay between convection and nucleosynthesis, so 
that variations in the C/0 ratio can be driven by pos- 
itive/negative changes in both carbon and oxygen abun- 
dances. For instance, while an increment of C/O is predicted 
due to the TDU, when HBB takes place the C/O ratio is 
expected either to decrease, as long as C is converted into N 
by the CN cycle, or even to increase if O is efficiently burnt 
in favour of N by the ON cycle (see Sect. 14.311 . 



^ See the WEB page www.iofle.rssi.ru/astro/conduct/ 

^ iESOPUS web- interface at|http://stev.oapd.inaf.it/aesopus| 
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Table 1. Input prescriptions of the TP-AGB models 



model 



low-T opacities: 
CNO variations 



mass loss 



symbol/lin(^ 



S06H 


no 


Straniero et al. (2006) 


A /dash 


S06C 


yes 


Straniero et al. (2006) 


o/long-dash 


BH 


no 


Blocker (1995) 


A/solid 


EC 


yes 


Blocker (1995) 


■/dot 



Symbols and line styles used in the plots 



temperature reached at the bottom of the external enve- 
lope, plus further information concerning the average con- 
tent of the ejecta, namely the helium mass fraction, and 
the C, N, O and Na enhancement/depletion factors, in 
terms of the quantities [X/Fe], where [X/Fe]=log(X/Fe)- 
log(X/Fe)Q. The last column shows the ratio between the 
average C-l-N-l-O abundance in the ejecta and the initial 
value, which is assumed to represent the chemical mixture 
at the epoch of the star's formation. 



Variations in CNO abundances affect the low- 
temperature opacities essentially in two ways, i.e. i) for 
T < 3 000 K by modifying the equilibrium molecular pat- 
tern depending on the C/0 ratio and, to a less extent, ii) 
for larger temperatures by changing the contributions of the 
CNO atoms to both the continuum and line opacity. A de- 
tailed discussion can be found in Marigo & Aringer (2009; 
see their section 4.2) and Ventura & Marigo (2009). 



3.4 Nuclear reactions 

The nucl ear network included i n the code is described in 
details in lVentura fc D'Antonal (|2005al ). The cross-sections 
of the 64 rea ctions considered ar e taken from the NACRE 
compilation (jAngulo et al.l Il999l). with the ex ception of 
^^N(p.'y)"0. taken from iFormicola et all |200i), and the 
three proton-capture reactions of the Ne-Na cycle, taken 
from lHale et al.l (12002!) for the '^^'!<ie{p,-yf^Na. reaction, and 
from iHale et al.l (|2004j" ) for the two p-captures by sodium 
nuclei. 



4 EVOLUTIONARY PROPERTIES 

4.1 An overview of the models 

The stellar models discussed here were followed from the 
pre-MS phase to almost the complete ejection of the exter- 
nal envelope. The initial chemical composition of the gas is 
assigned a total metallicity (mass fraction) Z = 0.001 and a 
degree of a— enhancement [of/F e] = -1-0.4, with the reference 
solar mixture taken from iGrevesse fc Sauvaj (119981 ). 

To investigate how much the results are affected by the 
the interplay between the use of the opacities accounting for 
the CNO variations and the mass-loss description we calcu- 
lated 4 sets of evolutionary models, designated with S06H, 
S06C, BH, and BC, which differ in the adopted prescrip- 
tions as outlined in Table [l] Specifically, we consider two 
formalisms for the mass loss, i.e. Straniero et al. (2006) and 
Blocker (1995), and two treatments of the low-T opacities, 
depending on whether the underlying chemical mixture is 
kept fixed or accounts for changes in the CNO abundances. 
In order to better disentangle the effects of each prescription, 
the models cover all four combinations of the two parame- 
ters, opacity and mass loss. 

The resulting physical and chemical properties of the 
TP-AGB models described above are presented in Tabled 
For each stellar mass we show the number of thermal pulses 
experienced by the star, the final core-mass, the maximum 



4.2 The interplay between mass loss and 
molecular opacities 

We can appreciate the qualitative effects of the different de- 
scriptions of mass loss and molecular opacities from Fig. (TJ 
where we show the evolution of a SMq model calculated 
according to the prescriptions listed in Table [T] to which 
we refer for the meaning of the various symbols. Core H- 
and He-burning phases are not included in this plot, that 
starts from the beginning of the TP-AGB phase. Each point 
marks the quiescent stage of pre-TP luminosity maximum. 
The solid lines are iso-mass locii, and connect the four evo- 
lutionary sequences at the stages when the total mass of the 
star has been reduced to 2.9, 2.7, 2.5 and 2Mq. 

In all four cases considered here the SM© models share 
a few common features, namely: i) they experience HBB, 
which is usually associated to an over luminositjij effect, and 
ii) they enter the domain of C-stars, as the surface C/0 
ratio increases above unity due to the TDU. At the same 
time significant differences arise. 

The S06H and S06C models experience a much weaker 
mass loss at the beginning, evolving at approximately con- 
stant mass for many TPs; this is at odds with the behaviour 
of BH and BC models, where an efficient mass loss deter- 
mines an earlier extinction of HBB and its overluminosity 
(following a rapid cooling of the envelope structure). This 
circumstance is seen in the maximum luminosity attained, 
which is ~ 0.2 dex fainter than in S06H and S06C sequences. 

The role played by the opacity treatment can be un- 
derstood by examining the evolution of models sharing the 
same description of mass loss, e.g. the S06H and S06C mod- 
els. In the latter, the rapid increase in the mass loss rate as 
soon as the surface C/O exceeds unity (clearly detectable 
as a jump in M) favours an earlier reduction of the mass 
of the external mantle, which again causes an earlier drop 
in the luminosity L (we have a ~ 0.1 dex difference in the 
maximum luminosity in this case). The drop in L is associ- 
ated to a lower temperature at the bottom of the convective 
envelope, i.e. less favourable conditions for HBB. As a con- 
sequence, we may conclude that, in general, using a mass 
loss description only mildly dependent on the luminosity, 
and/or neglecting the changes in the molecular chemistry 
in the opacity computations when C/0> 1, correspond to 
larger temperatures in the external mantle, in favour of a 
more efficient HBB. 



^ In quiescent stages TP-AGB models with HBB are brighter 
than expected by the classical core-mass luminosity relation 
(e.g.Boothroyd & Sackmann 1991) 
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Figure 2. Left panel: Evolution of the core mass as a function of the current stellar mass for models with initial masses 2.5, 3 and 3.5 
Mq . Each point on the tracks corresponds to quiescent CNO burning phase before the occurrence of each TP. Refer to Table [T] for the 
meaning of symbols. Right panel: The same as the left panel, but referring to the evolution of 3.5, 4.0 and 4.5 Mq models. 

Table 2. Relevant properties of AGB models 



M/M0 


NTP 


Mc/M@ 


logCT^rcD 


Y 


[I2c/Fe] 


[i^N/Fe] 


[i^O/Fe] 


[Na/Fe] 


R(CNO) 


BH models 


6.00 


24 


1.029 


8.05 


0.347 


-0.707 


1.325 


-0.369 


0.248 


0.929 


5.50 


28 


0.987 


8.05 


0.340 


-0.686 


1.350 


-0.468 


0.283 


0.946 


5.00 


33 


0.947 


8.02 


0.329 


-0.505 


1.458 


-0.410 


0.402 


1.205 


4.50 


34 


0.916 


8.00 


0.312 


-0.233 


1.669 


-0.112 


0.674 


2.012 


4.00 


33 


0.888 


7.97 


0.293 


0.090 


1.781 


0.191 


0.935 


2.850 


3.50 


28 


0.857 


7.94 


0.270 


0.259 


1.902 


0.464 


1.147 


4.074 


3.00 


24 


0.822 


7.88 


0.250 


0.796 


1.942 


0.678 


1.123 


5.518 


2.50 


22 


0.746 


7.49 


0.257 


1.805 


0.519 


0.985 


0.416 


11.204 


2.00 


20 


0.702 


7.05 


0.261 


1.564 


0.463 


0.622 


0.346 


6.015 


BC models 


6.00 


25 


1.028 


8.05 


0.348 


-0.706 


1.323 


-0.360 


0.252 


0.929 


5.50 


30 


0.983 


8.04 


0.335 


-0.667 


1.362 


-0.470 


0.285 


0.971 


5.00 


33 


0.946 


8.02 


0.329 


-0.527 


1.479 


-0.390 


0.418 


1.258 


4.50 


35 


0.916 


8.00 


0.311 


-0.175 


1.609 


-0.133 


0.614 


1.795 


4.00 


32 


0.886 


7.97 


0.292 


0.121 


1.758 


0.190 


0.912 


2.747 


3.50 


27 


0.856 


7.94 


0.266 


0.360 


1.872 


0.494 


1.104 


4.012 


3.00 


22 


0.821 


7.63 


0.250 


1.226 


1.336 


0.687 


0.397 


4.477 


2.50 


18 


0.742 


7.25 


0.257 


1.684 


0.516 


0.988 


0.396 


9.322 


2.00 


16 


0.677 


6.88 


0.265 


1.710 


0.481 


0.897 


0.354 


9.065 


1.50 


12 


0.648 


7.40 


0.251 


0.310 


0.044 


0.400 


0.030 


1.129 


S06H models 


6.00 


73 


1.044 


8.07 


0.355 


-0.320 


1.648 


-0.608 


0.014 


1.743 


4.50 


50 


0.924 


8.01 


0.320 


0.298 


2.023 


0.001 


0.648 


4.377 


4.00 


46 


0.900 


7.99 


0.304 


0.404 


2.133 


0.193 


0.948 


5.729 


3.50 


36 


0.883 


7.97 


0.270 


0.360 


2.216 


0.373 


1.184 


6.997 


3.00 


39 


0.841 


7.93 


0.261 


0.677 


2.457 


0.696 


1.837 


12.570 


2.50 


52 


0.806 


7.89 


0.270 


0.739 


2.677 


0.983 


2.436 


21.029 


2.00 


50 


0.784 


7.82 


0.271 


1.268 


2.652 


1.033 


1.940 


22.075 


1.50 


28 


0.731 


7.03 


0.260 


2.033 


0.156 


0.948 


0.385 


16.285 


S06C models 


5.00 


54 


0.956 


8.03 


0.334 


0.018 


1.856 


-0.250 


0.377 


2.907 


4.50 


47 


0.921 


8.01 


0.317 


0.083 


1.973 


-0.059 


0.585 


3.825 


4.00 


40 


0.890 


7.98 


0.295 


0.389 


2.109 


0.283 


0.997 


5.580 


3.50 


21 


0.846 


7.95 


0.270 


0.251 


2.202 


0.490 


1.310 


6.998 


3.00 


31 


0.830 


7.92 


0.251 


0.675 


2.270 


0.747 


1.518 


9.168 


2.50 


23 


0.758 


7.51 


0.253 


1.833 


0.510 


1.113 


0.428 


12.840 


2.00 


21 


0.702 


7.03 


0.260 


1.534 


0.476 


0.621 


0.338 


5.714 


1.70 


18 


0.685 


6.88 


0.253 


1.573 


0.053 


0.578 


0.059 


5.909 


1.50 


15 


0.675 


6.76 


0.253 


1.388 


0.107 


0.469 


0.082 


4.040 


1.20 


14 


0.660 


6.57 


0.253 


1.141 


0.056 


0.412 


0.056 


2.608 


1.00 


13 


0.671 


6.37 


0.259 


2.131 


0.073 


1.311 


0.144 


23.516 



4.3 C-star stage and HBB quenching 

Figure [5] shows the evolution of the C-O core mass during 
the TP-AGB phase of models with different initial masses 
and/or input proscriptions. The BH and BC models (full 
squares and open triangles) evolve to smaller core masses, 
as a result of the higher mass loss in their early TP-AGB 
phase. 

In the left panel we see that in the less massive BH 



and BC models the mass of the remnant is almost indepen- 
dent of the adopted opacity, being only slightly higher, a s 
expected, in the BH case. When the lStraniero et al.l (|2006r ) 
mass loss rate is used (S06H and S06C models) the masses of 
the remnant differ more significantly (see col. 3 of Table [5)), 
the discrepancy consisting in 5Mc ~ O.OSMq for M — 2Mq, 
SMc ~ O.O4M0 for M = 2.5Mq, and 5Mc - 0.02Mq for 
M — 3.5 Mq . The difference in Mc is higher at lower masses, 
because these models reach more easily the C-rich stage, 
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Figure 3. The same as Fig. [2] but siiowing ttie temperature at tlie bottom of the convective envelope 



that is accompanied by an increase in the molecular opacity 
in the S06C models. The stronger sensitivity of the S06H 
and S06C models to the opacity treatment compared to 
the BH and B C models can be explained by the fact the 
iBlockerl ([1995|) formalism predict quite high mass-loss rates 
so that the models become C-rich when most of their enve- 
lope masses have already been lost, with consequent little 
possibility for establishing great differences in the mass of 
the remnants. 

The evolution of Mc for higher mass models is shown 
in the right panel of Fig. (2] The BH and BC lines are again 
almost over-imposed in this diagram, and the differences 
between the S06H and S06C sets of models persist up to 
masses of the order of 4Mq , whereas in more massive models 
the surface C/0 ratio hardly approaches unity, thus render- 
ing the results almost opacity-independent. For the 4.5Mq 
model we note that the tracks are practically split into two 
branches, according to the mass loss treatment. 

The description given above can be complemented with 
the behaviour of the temperature at the bottom of the con- 
vective envelope, Tbco- This is the key-quantity to under- 
stand the degree of nucleosynthesis expected, and thus the 
chemical composition of the ejecta. The evolution of Tbco as 
a function of the current stellar mass is shown for all models 
in Fig. [3 

For M ^ 2.5Mq HBB conditions, i.e. Tbce > 60 x lO'' 
K, are never reached, with the only exception of the S06H 
models, where the rela tively mild mass loss, due both to the 
IStraniero et al.l (|2006l ) prescription itself and to the molec- 
ular opacities that neglect CNO variations, favours the in- 
crease in the core mass and temperature at the bottom of 
the envelope, eventually reaching HBB conditions. 

This is confirmed by the two panels of Fig. 2] show- 
ing the evolution of the surface carbon abundance for the 
two models with initial masses 2.0 (left panel) and 2.5 Mq 
(right panel). In models BH (solid line) and BC (dotted line) 
the surface carbon increases with a step-like profile that is 
typical of TDU effects. The BC models undergo a smaller 
number of TPs as a consequence of the higher mass loss 
rates experienced as the C/0 ratio exceeds unity, therefore 
the surface carbon is lower. 

An increasing trend of the carbon abundance is also 
shown by the S06C models (long-dashed), although in this 
case most of the carbon enhancement is achieved at the be- 
ginning, when the mass lost is negligible; at later stages the 
strong mass loss triggered by the formation of C-bearing 
molecules prevents further meaningful changes in the sur- 
face carbon mass fraction. 

Different is the behaviour of the S06H models (dashed 
lines) in both panels of Fig. 2] the fact that HBB is operating 
can be recognized from the rapid drop of the surface carbon 
abundance (being converted into nitrogen). In these cases 
the ejecta of the stars are expected to be carbon-depleted. 

Moving to higher masses, we see from Fig. 14.31 that in 
the 3Mq model HBB conditions are reached in both the 
S06C a nd S06H mode ls, whereas in the models calculated 
with the lBlocken (| 19951 ) mass loss HBB is quenched when the 



correct opacities are adopted (compare the dotted and solid 
tracks in Fig. 14.31 and the corresponding evolution of Tbcc 
marked by open triangles and full squares in the bottom- 
left panel of Fig. [3]). HBB is found in all models of higher 
mass, as can be seen in the right panels of Fig. [31 The only 
difference we see here is the different mass left in the en- 
velope when the asymptotic temperature at the bottom of 
the envelope is reached. For the same reasons previously dis- 
cussed, S06H and S06C models reach this stage when only 
a small amount of mass has been lost: this will have some 
effects on the yields (see Sect.[S]). 



4.4 Summary 

The results confirm the analysis bv lMarigol (1200211 a nd th e 
main findings of the investigation by ICristallo et al.l l|200a ) : 
low-mass AGB models need to be calculated with the cor- 
rect low-T opacities, that account for the changes in the 
molecular chemistry driven by variations of the C/0 ratio, 
and more generally in the CNO abundances. Neglecting the 
surface carbon enrichment in the opacity computations de- 
lays the ejection of the external mantle, leads to a higher 
final core mass and, in more massive models, leads to higher 
temperatures at the envelope base, hence favouring the de- 
velopment of HBB. 

The extent of the changes introduced by a correct opac- 
ity computation is conditioned to the treatment of mass loss. 
The indirect effects of using C-rich opacities on the growth 
of the remnant core, maximum luminosity and temperature 
at the bottom of the convective envelope are larger when the 
mass loss in the early AGB phase is sufficiently mild. In fact, 
under these circumstances, the C-rich stage is reached when 
the envelope still contains a large fraction o f mass; this ap- 
plies, for instance, when a treatment like the lStraniero et ahl 
(2006) is used. Models calculated with the iBlockej (Il995t ) 
formula are less sensitive to the opacity changes. 

The analysis on HBB made by IMarigo (2007 ') is con- 
firmed. For the reasons mentioned a bove, the quenc hing of 
HBB in models calculated with the iBlocken (|l995r ) mass- 
loss description is restricted to a narrow range of masses 
(clustering around 3M0 in the present investigation). Con- 
versely, when a more modera te mass loss, e.g . the V W93 or 
the modification suggested bv lStraniero et al.l l|2006l) is used, 
the extinction of HBB involves a wider range of masses, that 
in this investigation spans the interval from 1.8 to 2.6Mq. 

All models with M^ 3.5M0 achieve HBB conditions 
regardless of the opacity treatment, because the C/O ratio 
keeps in any case well below unity. Some differences between 
models differing in the mass loss treatment persist up to 
~ 4.5A/0, and tend t o disappear for larger masses. As ex- 
tensively discussed in I Ventura fc D'Antonal (|2005al ). in this 
range of masses the key quantity, that controls the activation 
and the strength of HBB, is the treatment of the convective 
instability, and, specifically, the efficiency of convection. 
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Figure 4. Left: the variation of tlie surface carbon abundance in models with initial mass 2M0, calculated for various treatments of 
mass loss and molecular opacities. Right: the same as in the left panel, but for models with M=2.5M0. The meaning of the different 
tracks is as follows, solid: BH models; dotted: BC; dashed: S06H; long-dashed: S06C 



Figure 5. The CNO variation in the ejecta for the 4 sets of 
models under consideration. The ordinate shows R(CNO), i.e. the 
ratio between the average C+N+O in the ejecta and the initial 
C+N+O content of the star. The meaning of the different lines is 
the same as in Fig. |4] 



5 HOW ROBUST ARE THE AGB YIELDS? 

Changes in the surface chemical composition of AGB stars 
are caused by HBB and TDU. The former keeps the overall 
C+N+O abundance constanlQ whereas TDU increases the 
carbon abundance, hence the total C+N+O, independently 
of the possible later conversion (via HBB) of ^^C to '^^N at 
the bottom of the convective zone. 



5.1 The C+N+O increase 

In the context of the distinctive chemical patterns exhib- 
ited by GC stars, the surface C+N+O is a key-quantity 
in the debate concer ning the self-enrichment scenario by 
mass ive AGB stars (JD'Ercole et al.l l2008l : IVentura et al.l 
[200l|), because the spectroscopic surveys of GC stars con- 
firmed that the t otal C+N+O is approximately constant 
Ijlvans et al.lll999l ). even when comparing stars with differ- 
ent surface abundances of other elemental species. There- 
fore, AGB ejecta with a great CNO enhancement would rule 
out their role as the possible cause of the observed patterns 
l|Karakas et al.ll2006l '). 

The CNO increase in the ejecta is shown in Fig. (5] this 
is the same quantity reported in col. 10 of Tab. (2] In all 
cases R(CNO) shows a ma^ximum for M~ 2 — 2.5Mq: lower 
masses experience fewer TDUs, whereas in more massive 
models the effects of TDU are reduced, due to dilution of 
the dredged-up, C-rich material into a massive envelope. 

For M > 3.5Mq, similarly to other quantities discussed 
in Sect. |4l R(CNO) depends mainly on the mass-loss pre- 
scription; the models split into two branches. S06H and S06C 
ejecta show a greater enhancement in the total CNO, due 
to the larger number of TPs and TDUs experienced. Inter- 
estingly, for M~ 6M0 the two branches tend to converge 
to unity (i.e. invariance of the total C+N+O), since these 
models are characterised by quite efficient HBB and large 
luminosities, which lead to a quick ejection of the envelope. 

For Ad < 3.5A'/q the treatment of the molecular opac- 
ities becomes more important: the CNO enhancement is 
lower for S06C and BC models compared to S06H and BH. 
The largest R(CNO) is reached by the S06H models (dashed 
line) due to the great number of TPs experienced, while the 
others (BH, BC and VW models) show more similar trends. 

The left panel of Fig. |B] shows the chemical content of 
the ejecta, on the [C/Fe] vs. [N/Fe] plane. Models with the 



^ We recall that the invariance applies to the total C+N+O abun- 
dance expressed as number fraction, and not as mass fraction: in 
fact, the CNO-cycle preserves the total number of the CNO cat- 
alysts, not their total mass. 



same mass are connected with continuous lines. The points 
in the lower-right corner correspond to low-mass models 
(M ^ 2.5Mq) not experiencing any HBB, and are char- 
acterized by large [C/Fe] as a consequence of the TDU. In 
general, a large value of [N/Fe] is a signature of HBB. The 
spread in [N/Fe] relative to models with masses in the range 
2-3M0 is a mere consequence of the different (if any) de- 
gree of HBB achieved in these masses. S06H models (full 
triangles) constitute a sort of upper envelope for the N en- 
richment: as discussed previously, these models experience 
a large number of TDUs, and HBB is always efficient. 

Note the position of the I.SMq BC model, indicated 
by the full square with the lowest [N/Fe] : the chemical con- 
tent of the mass ejected shows only a modest increase in 
carbon and no change in nitrogen, because the higher mass 
loss quickly leads to end of the AGB phase. This is also 
confirmed in Fig. [S] (right panel) by the position of the 
same point in the [0/Fe]-[Na/Fe] plane, corresponding to 
[O/Fe]=+0.4 and [Na/Fe]=0, i.e. the initial chemical com- 
position assumed for these stars. 

Higher mass models experience a more powerful HBB, 
thus exhibit a lower [C/Fe] in their ejecta. The relatively 
modest [N/Fe] is due to the lower efficiency of TDU in in- 
creasing the surface abundance of carbon when the dredged- 
up material is diluted into a massive envelope: for the same 
reason this limits the increase of nitrogen, which is produced 
by proton-captures on carbon nuclei. 



5.2 Sodium 

The right panel of Fig. [S] shows the oxygen and sodium con- 
tent of the yields, in terms of [Na/Fe] vs. [O/Fe]. These 
two elements are the main targets in the spectroscopic sur- 
veys of GCs, and the oxygen-sodium anti-correlation is a 
well-established pattern observed in practically all the GCs 
(|Carrettall2006l ). 

The ejecta of the lowest masses are 0-enriched (we re- 
call that the initial abundance of our a— enhanced mixture is 
[O/Fe]=+0.4), due to TDU. The S06C models (open circles) 
show a slightly higher [O/Fe], because they experience many 
TDUs, and HBB is far from being activated, thus keeping 
the surface oxygen high. The greater is the mass, the lower is 
the oxygen content of the ejecta, due to HBB which destroys 
^''O via proton captures in favour of '^4N. 

The behaviour of so dium at the surface of AGB stars 
is described in details in I Ventura fc DAntonal (|2008l ') (see 
their figure 4). The surface sodium content is unchanged in 
the lowest masses (see also col. 8 in Tab.m, while increases 
in models (M^ 2M0) where the temperature at the bottom 
of the envelope is high enough to convert the dredged-up 
^^Ne to Na . Similarly to nitrogen, the largest differences in 
[Na/Fe] due to mass-loss and opacity prescriptions are found 
in those models with mild HBB (note again the higher values 
attained by the S06H models). For higher temperatures the 
destruction channel of Na via proton capture prevails over 
production, leading to the lower [Na/Fe] for larger stellar 
masses. 
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Figure 6. Carbon and nitrogen (left panel), and oxygen and sodium (right panel) contents in the AGB ejecta. Lines connect models 
with the same initial mass. 



5.3 An overview on the uncertainties in the yields 
by AGB stars 

The analysis developed so far allows us to have a first ap- 
praisal of the uncertainties alTecting the yields from AGB 
stars produced by two main factors: mass loss and low-T ra- 
diative opacities. We have explored two treatments of mass 
loss which, especially in the early phases of the AGB evolu- 
tion, differ considerably: in this way our investigation should 
sample the present status of indetermination in the chemical 
yields. 

In agreement with previous investigations, we confirm 
that using opacities coupled to the actual surface composi- 
tion is mandatory for masses below a threshold value, corre- 
sponding to the full activation of HBB, which in the present 
work corresponds to M ~ 3.5 Mq. When the changes in 
the molecular chemistry are correctly accounted for in the 
opacities, the attainment of the C-rich stage (C/0> 1) is 
usually accompanied by a sudden increase in the mass loss, 
which favours the ejection of the whole mantle, speeding up 
the end of the AGB evolution. The enrichment in carbon 
and CNO is consequently lower. 

Typically, the changes in the yields, determined by use 
of the correct opacities, are more significant in models where 
a milder dependency of mass loss on luminosity is adopted. 
In fact, this would correspond to a longer TP-AGB lifetime 
during which the surface composition may be affected by i) 
a larger number of TDU episodes, and ii) more favourable 
conditions for the development of HBB. Conversely, assum- 
ing a ve ry efficient mass-loss formalism, as predicted by a 
iBlocken (l995)-like law, the accuracy loss introduced by un- 
suitable opacities (e.g. valid for scaled-solar mixture) is less 
dramatic, since the carbon-rich stage is reached when a con- 
siderable fraction of the mass of the envelope is already lost. 

Within this general trend, the sensitiveness of the yields 
to the model assumptions vary according to the mass of the 
star, as outlined below. 

The largest uncertainties in the yields belong to stars 
for which the temperature at the bottom of the convective 
envelope is just sufficient to ignite HBB, i.e. 1.8 — 3.0 Mq. 
In particular, [N/Fe] and [Na/Fe] may vary by '^ 2 dex, 
depending on the assumed input prescriptions. The oxygen 
yield is more robust. 

At lower masses, i.e. M ~ l.SM©, the results are more 
affected by the mass -loss descriptio n: the TDU does not even 
take place when the lBlockeij l| 19951 ) prescription is followed, 
being prevented by a rapid expulsion of the external enve- 
lope. 

Above 3.5 Mq the opacity treatment become less influ- 
ential, and most of the uncertainty in the yields should be 
ascribed to the the mass-loss treatment. 

This investigation has also shown that, although the 
yields may vary according to the adopted mass-loss rate in 
any range of mass, the predicted trends between different 
elemental species do not depend on the model assumptions, 
e.g. the slopes in both [C/Fe]-[N/Fe] and [O/Fe]- [Na/Fe] di- 
agrams are preserved (see Fig. [6]). 



Figure 7. Mass-loss rates experienced by a 2Mq model, com- 
puted with the ATON code, according to various mass-loss pre- 
scriptions available in the literature. 



6 A COMPARISON WITH PREVIOUS 
INVESTIGATIONS 

In view of testing the reliability of AGB models, we compare 
our results with investigations by other r esearch groups, fo- 
cusing on t he the recent investiga tions bv lWeiss &: FergusonI 
(2009) and ICristallo et al.l (|2009l ). in which low-T opacities 
are suitably calculated by considering the increase in the 
surface carbon abundance, an essential requirement to cor- 
rectly describe C-rich stage. 

We consider, in particular, the models computed with 
same initial mass, 2.0 Mq, that include our S06C and 
BC models, the Z = 0.001 model by ICristallo erall 
(.2009i') (see their Tab . 3), an d the Z = 5 x 10"" model 
bv IWeiss fc FergusonI (|2009l ) (see their Tab. B.4). Al- 
thou gh the init ia l chem ical compositions are not the same, 
( Cris tallo et al.l 1120091) use a scaled-so l ar mi xture with 
Z = 0.001, whereas iweiss fc FergusonI (12009!) adopt an 



Q— enhanced mixture with Z=0.0005), the main differences 
can be explained in terms of the different input prescrip- 
tions. 

One striking featur e is the different num ber of TPs ex- 
perienced, that is 19 in lCristallo et al.l (|2009|). 21 and 16 in 
our tw o S06C and BC cases, and only 6 in lWeiss fc FergusonI 
(J200g|) , despite of the fact that this latter model corresponds 
to the lowest metallicity. The small number of TPs found 
in the Weiss & Fergusonj (|2009|) work is accompanied by a 
shorter AGB life-time, which is ~ 1/3 of our estimated val- 
ues. Due to the absence of HBB, we rule out convection as a 
possible cause of such a discrepancy, and focus on the other 
two, uncertain, physical inputs, i.e. mass loss and the treat- 
ment of the convective borders, particularly during the TDU 
phase. 

Let us first focus on the effect produced by differ- 
ent mass-loss treatments. Figure. [7] shows the evolution of 
the mass-loss rate in our 2.0 Mq AGB model computed 
with the ATON code, assuming CNO-var iable molecular 
opacit ies and three c hoices for M, n amely: IStraniero et al.l 
(120061) (open circl es), Blockerl ([1995J) (open triangles), and 
iWeiss fc FergusonI (|2009l ) (full squares). We recall that 
this latter scheme assumes a Reimers' formula up to a 
threshold pulsationa l period of 400 days, followed by the 
IWachter et al.l (|2008l ) rat e when the C-st a r stag e is reached. 

We note that the IStranie ro et al.l (|2006l ) treatment 
favours rates significantly smaller, particularly at the be- 
ginning of the AGB phase, with differences up to ~ 1.5 dex 
co mpared to the o thers. A smaller difference exists between 
the lBlockeij (| 19951 ) and the lWeiss fc FergusonI (|2009l ) recipes, 
the latter being higher by ~ 0.2 — 0.3 dex during most of 
the evolution. 

In the early AGB the IWeiss fc FergusonI (|2009l ) treat- 
ment predicts rates of the order of 10~^ — 10~® Mq/jt, which 
are, however, not sufficient to diminish the total number 
of TPs so as to recover their original result (6). This sug- 
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gests t hat the main reasons for the smaUer number of TPs 
in the IWeiss fc FergusonI l|2009f ) investigation compared to 
our models is an earlier transition to the C-star stage. This 
difference is due to two reasons: a) The modelling of con- 
vective overshoot from the bottom of the envelope (see their 
equation 2), based on an exponential decay of the diffu- 
sive coefficient, with the same e- folding distance used for 
main-sequences width fitting. This approach determines a 
much deeper extension of TDU when compared with our 
models (where no ove rshoot is assumed) , and with those by 
I Crist alio et al.l (|2009l ). who adopt an exponential decay of 
the convective velocity from the inner border of the external 
mantle, with an e-folding distance tuned in order to allow 
the synthesis of a given amount of ^''C available for neutron- 
captures. This amount of extra-mixing is also the plausible 
reason for the smalle r number of TPs experienced by the 
ICristallo et al.l (|2009f ) model compared to our S06C, and to 
the different final core-masses (0.67 Mq, to be compared to 
0.70 Mq); b) the convective overshoo t from the He-burning 
shell, adopt ed by Weiss fc FergusonI (|2009|), and neglected 
here and in [C ristallo etajj ( 20091 '). which, as described by 
iHerwig fc Aus tin (200j), renders the pulses stronger. The 
models presented here, though in cluding no over s hoot, are 
much more similar to those by ICristallo et al.l (|2009l ) in 
terms of duration of the whole AGB phase and of the chem- 
istry of the ejecta. 

A detailed investigation of the TDU phenomenology, 
and the possible effects of assuming some convective over- 
shoot from the bottom of the convective envelope, is be- 
yond the sc ope of the pre s ent paper; the interested reader 
ma y find in iHerwi g (20001. iHer^ (|2005l ), lMow"iwil (|l999l ) 
and lStraniero et al. (f997 ) a full discussion on the key issues 
relevant to address this problem. 

Anyhow, the present study clearly illustrates that the 
AGB evolution is the product of a complex interplay be- 
tween various processes: the treatment of the extra-mixing, 
for instance, has a remarkable impact not only on the chem- 
ical but also on the physical properties of these stars. 



7 CONCLUSIONS 

In view of improving the robustness of the results provided 
by AGB modelling (or at least to quantify the uncertainty 
range), we have investigated the sensitivity of the main phys- 
ical and chemical properties to the choices inherent the input 
macro- and micro-physics. 

We find a threshold mass, that is M ~ 3.5 Mq for 
Z — 0.001 in the present analysis, (note that it also depends 
on the assumed overshooting from the central regions dur- 
ing the core H- and He-burning), above which the results 
are little sensitive to the opacity treatment, because HBB 
prevents the surface C/O ratio to exceed unity. 

In this range of masses mass loss modelling plays a 
crucial role. Models calculated with a treatment of the 
mass loss only mildly dependent on the lu minosity, such as 
the cla ssic I Vassiliadis fc WoodI (119931 ) or the lStraniero et all 
(|2006r) recipe used here, predict smaller rates during most 
of the AGB phase, thus more TPs are experienced. This 
also favours the growth of the luminosity and temperature 
at the bottom of the convective zone, and a more advanced 
HBB nucleosynthesis: the yields will contain the signature of 



proton-capture nucleosynthesis, e.g. lower carbon and oxy- 
gen abundances, and higher nitrogen content. However, the 
general patterns of the C-N and 0-Na relations, and the very 
small increase in the overall CNO abundance, are almost in- 
dependent of the opacity and mass-loss treatment: we expect 
convection modellin g plays a relevant role here, in a greement 
with the analysis bv lVentura fc D'Antonal (|2005al ). 

For masses smaller than 3.5 Mq the interplay between 
opacity and mass-loss treatment is more tricky, because the 
lack of HBB favours the increase in s urface carbon. Gener- 
ally speaking, we confirm th e resu lts by Marigol ( 20021 . : 



2007), 



20091 ) 



ICristallo etal] (|2008| . I2009D and IWeiss fc FergusonI 
the use of the correct opacities in the low-T regime speeds- 
up the loss of the stellar envelope, via the strong winds as- 
sociated to the general expansion of the structure when the 
surface C/O ratio exceeds unity. We expect in this case a 
smaller number of TPs, and a less advanced nucleosynthesis 
(if any) at the bottom of the convective envelope. 



Models calculated with a mass-loss rate steeply depen- 
dent on the luminosity show a more homogeneous behaviour, 
whereas when a milder dependency of M with luminosity, 
e.g. a relationship that relates M to the stellar period, is 
adopted, many physical properties, e.g. core-mass, temper- 
ature at the bottom of the convective zone, total number 
of thermal pulses experienced, follow a completely different 
behaviour. In this case the quenching of HB B determ i ned b y 
the use of the C-rich opacities, predicted bv iMarigol (|2007|), 
is confirmed, being a common feature of all the models with 
1.8 — 3.0 Mq. a narrower range of masses (around ~ 3 Mq) 
is otherwise affected by this un certainty when an efficient 
mass-loss prescription, like the IBlockeil ([1995|) formula, is 
adopted. 



Contrary to the models experiencing HBB, the yields 
from lower mass stars are highly uncertain and model- 
dependent: the average [C/Fe] in the ejecta shows up an 
overall uncertainty of 0.7 dex, though this difference reduces 
to ~ 0.3 — 0.4 dex if models calculated with the same opac- 
ities are compared. The situation for nitrogen and sodium 
is more extreme, as they are particularly sensitive to the 
ignition of HBB. The differences in [N/Fe] and [Na/Fe] can 
reach ~ 2 dex, and even when comparing models with the 
same assumptions for the opacity treatment, discrepancies 
of one order of magnitude still persist. The oxygen yields 
appear, instead, more stable. 

These trends should be considered representative of 
AGB models with a low metal content {Z ~ 0.001) while 
they need to be further explored at larger metallicities i.e. 
including those characteristic of the LMC (Ventura et al., in 
preparation). Moreover, in the present unresolved scenario 
about the efficiency of TDU and mass loss, a valuable con- 
tribution may come from synthetic AGB models, i.e. via a 
thorough calibration of AGB properties as a function of stel- 
lar mass and metallicity based on a large set of observables 
(HR diagrams including the MID infrared, C/M star ratios, 
chemistry of atmospheres and planetary nebulae, initial-final 
mass relation, integrated colors and brightness fluctuations). 
This work is in progress. 
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